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Abstract 

Here presented is a unified approach to Stirling numbers and their general- 
izations as well as generalized Stirling functions by using generalized factorial 
functions, fc-Gamma functions, and generalized divided difference. Previous 
well-known extensions of Stirling numbers due to Riordan, Carlitz, Howard, 
Charalambides-Koutras, Gould-Hopper, Hsu-Shiue, Tsylova Todorov, Ahuja- 
Enneking, and Stirling functions introduced by Butzer and Hauss, Butzer, 
Kilbas, and Trujilloet and others are included as particular cases of our gen- 
eralization. Some basic properties related to our general pattern such as their 
recursive relations and generating functions are discussed. Three algorithms 
for calculating the Stirling numbers based on our generalization are also given, 
which include a comprehensive algorithm using the characterization of Riordan 
arrays. 
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1 Introduction 

The classical Stirling numbers of the first kind and the second kind, denoted by 
s(n, k) and S{n, k), respectively, can be defined via a pair of inverse relations 

n n 

Nn = I]s(n,fc).^^ z" = ^5(n,fc)[z]fe, (1.1) 

fc=0 fc=0 

*This paper was presented in the International Conference on Asymptotics and Special Func- 
tions, City University of Hong Kong, 30 May - 03 June, 2011 
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with the convention s(n, 0) = 5(n, 0) = dnfi, the Kroneckcr symbol, where z G C, 
n e No = N U {0}, and the fahing factorial polynomials {z]n = z{z~l) ■ ■ ■ {z — n + 1). 
\s{n,k)\ presents the number of permutations of n elements with k disjoint cycles 
while S{n, k) gives the number of ways to partition n elements into k nonempty sub- 
sets. The simplest way to compute s{n, k) is finding the coeflacients of the expansion 
of [z]n- [20] gives a simple way to evaluate S{n, k) using Horner's method. 

Another way of introducing classical Stirling numbers is via their exponential 
generating functions 

n>/c n>k 

where < 1 and k G Nq. In [26], Jordan said that, "Stirling's numbers are of the 
greatest utility. This however has not been fully recognized." He also thinks that, 
"Stirling's numbers are as important or even more so than Bernoulli's numbers." 

Besides the above two expressions, the Stirling numbers of the second kind has 
the following third definition (see [13] and [26]), which is equivalent to the above two 
definitions but makes a more important rule in computation and generalization. 



= ^E(-iv(-)(fc-^r. (1.3) 

Expressions (1.1) - (1.3) will be our starting points to extend the classical Stirling 
number pair and the Stirling numbers. 

Denote {z)n,a '■= z{z + a) ■ • ■ {z + {n— l)a) for n = 1, 2, . . ., and {z)o^a = 1) where 
{z)n,a is called the generalized factorial of z with increment a. Thus, {z)n,-i = [z]n 
is the classical falling factorial with [^;]o = 1, and {z)n,o = z". More properties of 
{z)n,a will be presented below. 

With a closed observation, Stirling numbers of two kinds defined in (1.1) can be 
written as a unified Newton form: 

n 
k=0 

with S{n, k, 1, 0) = s{n, fc), the Stirling numbers of the first kind and S{n, k, 0, 1) = 
S{n,k). the Stirling numbers of the second kind. Inspired by (1.4) and many exten- 
sions of classical Stirling numbers or Stirling number pairs introduced by [8], [23], 
[46], [24], etc. We may define a unified generalized Stirling numbers S{n,k,a,f},r) 
as follows. 

Definition 1.1 Let n G N and a,/3,r G M. A generalized Stirling number denoted 
by S{n, k, a, /3, r) is defined by 
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{z)n-a = Y^S{n,k,a,p,r){z -r)k-i3. (1.5) 

fe=0 

In particular, if {a,0,r) = (1,0,0), S'(n, fc, 1, 0, 0) is reduced to the unified form of 
Classical Stirling numhers defined by (1-4-) ■ 

Each {z)n,-a does have exactly one such expansion (1.5) for any given z. Since 
deg {z — r)k-i3 = k for all k, which generates a graded basis for 11 c F — )• F, the 
hncar spaces of polynomials in one real (when F = M) or complex (when F = C), in 
the sense that, for each n, {{z — r)„,_^} is a basis for n„ c 11, the subspace of all 
polynomials of degree < n. In other wards, the column map 

W,:¥^ ^U:s^Yl ^("' ^' 0(^)^.-/3' 

fe>0 

from the space Fq^ of scalar sequences with finitely many nonzero entries to the 

space n is one-to-one and onto, hence invertible. In particular, for each n G N, the 
coefficient c(n) in the Newton form (1.5) for {z)n-a depends linearly on {z)n-a-, 
i.e., {z)n-a ^ s{n) = (W~^(2;)„,_„)(n), the set of S{n,k,a, j3,r), is a well-defined 
linear functional on 11, and vanishes on n<„_i. 

Similarly to (1.1), from Definition 1.1 a Stirling-type pair {3^,8"^} = {S^(n,k), 
S'^{n,k)} = {S{n,k; a,l3,r), S{n,k;P,a,—r)} (see also in [24]) can be defined by 
the inverse relations 



{z)n-a = ^S'^{n,k){z-r)k,-0 

n 

{z)r,,-p = J2s^n,k){z + r)k,-a, (1.6) 

fc=0 

where n G N and the parameter triple (a, /3, r) ^ (0, 0, 0) is in or C'^. Hence, we 
may call S^ and S"^ an (a, /3, r) and a {13, a,— r)— pair. Obviously, 

^(n. A:; 0,0,1) = 

because z" — X]fc=o (fc)('^ ~ l)*^- addition, the classical Stirling number pair 
{s{n,k), S{n,k)} is the (1,0,0)— pair {S^,S^}, namely, 

s{n,k) = S\n,k;l,0,0) S{n,k) = ^^(n, fc; 1, 0,0). 

For brevity, we will use S{n, k) to denote S{n, k, a, /3, r) if there is no need to indicate 
a, /?, and r explicitly. From (1.5), one may find 

5(0,0) = 1, 5'(n,n) = l, S(l,0) = r, and 5(n, 0) = (r)„,_„. (1.7) 

Evidently, substituting n — k — into (1.5) yields the first formula of (1.7). Com- 
paring the coefficients of the highest power terms on the both sides of (1.5), we 
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obtain the second formula of (1.7). Let n = 1 in (1.5) and noting S{1, 1) = 1, we 
have the third formula. Finally, substituting z = r in (1.5), one can establish the 
last formula of (1.7). The numbers a(n,k) discussed by Doubilct et al. in [16] and 
by Wagner in [47] is k\S{n,k;0, 1,0). More special cases of the generalized Stirling 
numbers and Stirling- type pairs defined by (1.5) or (1.6) are surveyed below in Table 
1. 



(a,/3,r) 


dual of {a, j3, r) 


S{n, k) pairs 


Name o f Stirling numbers 


(-1.1.0) 


(1.-1.01 


{-l)"-^'n\{'r])/kl 


L(iJ) iiiiirihrr pair^2^^^ 


(-1,0,0) 


(0,-1,0) 


\s{n,k)\ 
i-XY^-^Sin k) 


signless Stirling numbers[37] 


{1,0,O){0^O) 


(0,1,0) 


5(n,fc, 1,61,0) 
S{n,k,9, 1,0) 


Carlitz's degenerate Stirling 

number pair\7'\ 


(1,0, -A) 


(0,1, A) 


5(n,fc, 1,0,-A) 


Carlitz's weighted Stirling 

iiLifftuaf puii lOj 


(1,^,-A) 


(0,1, A) 


S'(n,fc, 1,6',-A) 

Ol/i, A/, C, ±, /\J 


Howard! s weighted degenerate 

kJt/ll Itliy llUblliUei pU/ti [■^'Jj 


(0,1, -a + 6) 


(1,0, -6+a) 


S'(n,fc, 0,1, -a + 6) 
o (^n, K, i,U, —0 + aj 


Gould — Hopper's non — central Lah 

number pair [18] 


(l/s,l,-a + 6) 


(l,l/s,-6 + a) 


^(n, fc, 1/s, 1, —a + b) 
^(n, fc, 1, 1/s, — 6 + a) 


Charalambides — Koutras's non— 
central C number pair [9, 10] 


(1,0, 6 -a) 


(0,l,a-&) 


S{n, k, 1, 0, 6 — a) 
S{n, k, 0, 1, a — b) 


Riordan's non — central Stirling 

number pair [34] 


{a, 13,0) 


(^,a,0) 


Aap{r,m) 


Tsylova's Stirling number pair[46] 


(a,/3,r) 


{P,a,-r) 


S{n, k, a, l3, r) 
S{n, k, l3, a, — r) 


Hsu — Shiue's Stirling 
number pair [24] 


(l,a:,0) 




ank{x) 


Todorov's Stirling number s[4:5\ 


(-l/r,l,0) 




B{n,r, k) 


Ahuja — Enneking's associated 
Lah number s[31\ 


(-1,0, r) 




S{n — r,k — r, —1, 0, r) 


Br Oder's r — Stirling number s[3] 







Table 1. Some generalized Stirling Numbers and Stirling Number pairs 



The classical falling factorial polynomials [z]„ = z{z — 1) • ■ • (z — n + 1) and 
classical rising factorial polynomials [z]"" = z{z + 1) • • • (^ + n — 1), 2: € C and n € N, 
can be unified to the expression 

{z)n,±i ■.= z{z±l)---{z±{n-l)), 
using the generalized factorial polynomial expression 
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{z)n,k-=z{z + k)---{z+{n-l)k) = {z + {n-l)k}n,-k {zeC,neN). (1.8) 

Thus {z)n.i = [z]" and = [z]n- In addition, we immediately have the rela- 

tionship between [z]" and {z)n,k as 

= {zGC,n€n,k>0). (1.9) 

Similarly, we obtain 

{z)„,-k = z{z-k)---{z-{n-l)k) = k''[z/k]n (^ G C,n e N, A; > 0). (1.10) 

The history as well as some important basic results of the generalized factorials can 
be found in Chapter II of [26], and an application of the generaUzed factorials in the 
Lagrange interpolation is shown on Page 31 of [17]. 

It is known that the falling factorial polynomials and rising factorial polynomials 
can be presented in terms of Gamma functions: [z]„ — T{z + l)/r(z — n + 1) and 
[z]"' = r{z + n)/T{z), and the gamma function r{z) can be defined in terms of 
factorial functions by (see, for example, [29]) 

r{z) = hm i^fi^ (^ e c - fcz_). (1.11) 

n^oo [Z\ 

As an analogy, the k-gamma function Tk, a one parameter deformation of the clas- 
sical gamma function, is defined by (see, for example [15]) 

Tkiz) := lim "■'^ (yfc>0,zeC-fcZ_). (1.12) 

n^oo {Z)n,k 

[z]"" and {z)n,k {k > 0) are also called the Pochhammer symbol and k-Pochhammer 
symbol, respectively. Even the parameter k is replaced by other parameters, we still 
call the corresponding Pochhammer symbol the /c-Pochhammer. 
For A; > 0, from (1.9), (1.11) and (1.12) (see also [27]) we have 

Tk{z) = fc(^/'=)-ir (I) . (1.13) 

Since [2;]" = r{z + n)/r{z), [4] extends the classical raising and falling factorial 
polynomials to generalized raising and falling functions associated with real number 
7 by setting 

■ r(z) ^^J^" r(2-7 + i)' ^ > 

respectively. We now extend {z)n,k defined by (1.8) to a generalized form associated 
with 7 e C using the relationship (1.9), namely, 

{z)^,k = k''[z/kr, (z)j,-k = k-^[z/k]^ (2GC,7eC,/e>0), (1.15) 
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which are caUed the generahzed raising and falling factorial functions associated 
with complex number 7, respectively. Using (1.13)-(1.15), we establish the following 
result. 

Theorem 1.2 If k > and {z)j^k is defined by (1.15), then 

There hold the following recurrence relations of the generalized raising and falling 
factorial functions. 

Proposition 1.3 If k > and {z)^^k is defined by (1.15), then there hold 



{z)^,k = {z + {l-l)k){z)^-i,k, {z)^-k = {z-{'y-l)k){z)^-i,-k. (1-17) 

In next section, we will present the unified expression and some properties of the 
generalized Stirling numbers of integer orders. Two algorithms based on the unified 
expression will be given. Then, we use the A:-Pochhammer symbol and A;-Gamma 
functions to extend the classical Stirling numbers of integer orders to the complex 
number orders in Section 3, which arc called the generalized, Stirling functions. The 
convergence and the recurrence relation of the generalized Stirling functions as well 
as their generating functions will also be presented. Finally, in Section 4 we will 
give more properties of generalized Stirling numbers and functions using the gen- 
erating functions of generalized Stirling numbers shown in Section 3, which include 
the asymptotic expansions of generalized Stirling numbers and functions and the se- 
quence characterizations of the Riordan arrays of generalized Stirling numbers. The 
third algorithm of the computation of the generalized Stirling numbers, including the 
classical Stirling numbers as a special case, will be shown using the characterizations 
of their Riordan arrays. 



2 Expressions of generalized Stirling numbers 

First, wc give an equivalent form of the generalized Stirling numbers S{n, k) defined 
by (1.5) by using the generalized difference operator in terms of /3 {j3 ^ 0) defined 

by 



A^/ = A^(A^-7) (fc>2) and A^/(t) /(t + /3) - /(t). (2.1) 

It can be seen that Ao(z)j^_^ = P'^k\6k.j, where 6h,j is the Kroneckcr delta 
symbol; i.e., Skj = 1 when k = j and otherwise. Evidently, from (1.10) there holds 



't' 




J. 


j 



I3''k\5k„ 



(2.2) 



Denote the divided difference of f(t) at t+i, i ~ 0,1, . . . ,k,hy f[t, t + 1, . . . , t+k], 
or [t,t + l, . . . ,t + k]f{t). Using the well-known forward difference formula, it is easy 
to check that 
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and 



^A''f{t) = f[t,t + l,...,t + k] = [t,t + l,...,t + k]f{t) 



^ A''0f{t) = f[t,t + p,t + 2p,...,t + kp] = [t,t + p,...,t + kp]f{t). 



pkkl 

We now give the following definition of the generalized divided differences. 
Definition 2.1 We define A^/(t) by 



BD\f{t) if 13 = 



where A^/(t) is shown in (2.1), f[t,t + p,...,t + kp] = [t,t + (3, . . . ,t + kp]f is the 
kth divided difference of f in terms of {t,t + 13, . . . ,t + and D''f(t) is the kth 
derivative of f{t). 

From the well-known formula 

D^fiO 



f[t,t + l3,t + 2p,...,t + kl3] 
where ^ is between t and t + k(3, it is clear that 



k\ 



which shows the generalized divided difference is well defined. 

Wc now give a unified expression of the generalized Stirling numbers in terms of 
the the generalized divided differences. 

Theorem 2.2 Let n, A; € No and the parameter triple (a, /3, r) ^ (0, 0, 0) is in 
or€?. For the generalized Stirling numbers defined by (1.5), there holds 



S{n,k,a,p,r)= Al{z)n,^ 



I z=r 



[r,r + p,...,r + k0\{z)n,-a ifl3¥^0 

if 13 = 0. 



(2.5) 



In particular, for the generalized Stirling number pair defined by (1.6), we have the 
expressions 
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S\n,k) = S\n,k,a,f3,r) = Aj(2)n,- 

[r,r + p,...,r + kp]{z)n,-a, ifPj^O 



S\n,k) = S^{n,k,p,a,-r) = A^(^)„,_^ 



-a ) 



if 13 = 



(2.6) 



j ^a{z)n-f3\^^_^ = [-r,-r + a,...,-r + ka]{z)n-fi, if a 

(2.7) 

Furthermore, if {a,0,r) = (1,0,0), then (2.5) is reduced to the classical Stirling 
numbers of the first kind defined by (1.1) with the expression 



s{n, k) = S{n, k, 1, 0, 0) = - ^'[^]„Uo • 

If{a,p,r) = (0,1,0), then (2.5) is reduced to the classical Stirling numbers of the 
second kind shown in (1.3) with the following divided difference expression form: 



S{n, k) = S{n, k, 0, 1, 0) = [0, 1, 2, ... , k]z\^^ . (2.8) 

The following corollary is obvious due to the expansion formula of the divided 
differences generated from their definition. 

CoroUciry 2.3 Let n, fc G No and the parameter triple (a, /3, r) ^ (0, 0, 0) is in 
or C^. If (3 ^0, for the generalized Stirling numbers defined by (1.5), there holds 



5(n,fc) = 5(n,fc,a,^,r) = -^XJ(-l)^Q(r + (fc-i)^)„,_a {n^O), (2.9) 
and 5(0, k) = dok- 

Remark 2.1 It can be seen from (2.9) that 

S{n, 0) = S{n, 0, a, /3, r) = (r)„,_„, (2.10) 

which is independent of (3 and has been shown in (1.7). The difference is deriving 
(2.10) from (2.9) needs {a,r) ^ (0,0) when ^ = 0. However, we have seen from 
(1.7) that the condition is not necessary. Another way to derive (2.10) using the 
characterization of the Riordan arrays of the generalized Stirling numbers will be 
presented in the Algorithm 4.5 in Section 4. 

Remark 2.2 If a/3 / 0, by taking the nth forward differences in terms of a and (3 
on the both sides of two equations of (1.6), respectively, one may obtain identities 
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n! 



n 
k=0 

The above two identities can be unified to be one: 

n 

n\a^ = J2 S{n, k, a, /3, r)A^{z - r)k,-0\^^^ . 
When a = 0, the above identity turns to 

n 

n\ = J2 S{n, k, 0, ^, r) D^{z - r)k,-pl^Q • 

fe=0 

Remarker 2.3 There exists another expression of the divided difference {z)n,-a 

in terms of Peano kernel of B-spline. Assume that the set r := {t,t + (3, . . . ,t + k^} 
lies in the interval [a, b] . Then on the interval, we have Taylor's identity 

\i (3 > 0, then is a weighted sum of values of derivatives of order < k, hence 
commutes with the integral in the above Taylor's expansion, which annihilates any 
polynomial of degree < k. Therefore, 



where 



M{y\T) := k[r, r + /3, . . . , r + - yf-^ 

is the Curry-Schoenberg B-spline (see [14]) with the knot set r and normalized to 
have integral 1. In particular, 

rb 

S{n,n,a,p,r) = A'l{z)n,-a = M{y\r,r + p, . . . ,r + np)dy = 1. 

We now present two algorithms for calculating generalized Stirling numbers. If 
/3 ^ 0, we denote 

^^f{t + £(3) := f[t, t + £l3,t+{£ + 1)13, ...,t+ jl3] (2.11) 

Thus, from (2.5) in Theorem 2.2, based on the recursive definition of tlic divided 
difference with respect to /3 (see Definition 2.1) 
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(2.12) 



we obtain an algorithm shown below. 



Algorithm 2.4 This algorithm of evaluating the generalized Stirling numhers is 
based on the construction of the following lower triangle array by using (2.11) and 
(2.12). 



{^)n,-<x\z=r 

{Z + I3)n-a\^^^ A^(z)„_„ 
{Z + 2/3) 

n, — a n, — a 



{z + kP)n,-al^, A^(z + (/c-l)/3)„,_„ Aj(z + (fc-2)/3)„,_ 

Table 2. The generalized Stirling numbers 



■ A. 



Thus, the diagonal of the above lower triangle array gives S{n,i,a, I3,r) = A^(^)„^_ 
for i = 0,1, . . . ,k. 

Example 2.1 We now use Algorithm 2.4 shown in Table 2 to evaluate the classical 
Stirling numbers of the second kind S{4:,k) = 5(4, /c, 0, 1, 0) {k = 1,2,3,4), which 
are re-expressed by (2.8). Thus, 












1 




1 




24 


= 16 


15 


7 


34 


= 81 


65 


25 


44 


= 256 


175 


55 



6 

10 1 

From the diagonal of the above lower triangular matrix, we may read 5'(4,0) = 0, 
5(4,1) = 1, 5(4,2) = 7, 5(4,3) = 6, and 5(4,4) = 1. Meanwhile, the subdiagonal 
gives 5(5, 1) = 1, 5(5, 2) = 15, 5(5, 3) = 25, and 5(5, 4) = 10. 

Example 2.2 For the Howard's weighted degenerate Stirling numbers 5(4, k) = 
5(4, fc, 1,1,-1), a similar argument of Example 2.1 yields 

(^)4,-iL_i=24 

(^ + l)4,-iL=_i =0 -24 

(^ + 2)4,-it,^_i =0 12 

{z + S)4,^i\l^_-^= -4 

(2 + 4)4,-1 =0 



Thus, 5(4, 0) = 24, 5(4, 1) 



1 

24, 5(4,2) = 12, 5(4,3) = -4, and 5(4,4) = 1. 
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Example 2.3 For the Howard's weighted degenerate Stirling numbers S'(4, k) = 
S(A, fc, 1, 2,-1), using Algorithm 2.4, we obtain S'(4, 0) = 24, 5(4, 1) = -12, 5(4, 2) = 
3, 5(4, 3) = 2, and 5(4, 4) = 1 reading from the following table. 

(^)4,-i|,=_i=24 

(z + 2}4,_i|,_i =0 -12 

(^ + 4)4,-i|,^_i =0 3 

(^ + 6)4-i|^^_i = 120 60 15 2 

+ 8)4^-1 = 840 360 75 10 1 

Remark 2.4 Obviously, Algorithm 2.4 is not limited to the case of /3 since 
when /3 = 0, A^(z)„ {k = 0, 1, . . . , n) on the diagonal of the lower triangle 

matrix in Table 1 are simply the l/k\ multiply of the derivatives D'^ {z)n,-a\ _ (see 
Theorem 2.2). 

Another algorithm based on the Horner's method can be established using a 
modified argument in the computation of the classical Stirling numbers of the second 
kind shown in [20]. More precisely, we have the following algorithm. 

Algorithm 2.5 First, we may write the generalized Stirling numbers 5(n, k) = 
S{n,k,a,p,r) defined by (1.5) (see Definition 1.1) as 



{z)n,-a = ^ 5(n, k){z - r)k-i3 



fc=0 

= 5(n, Q) + {z- r){S{n, 1) + (z - r - ^)(5(n, 2) + (z - r - 2/3)(5(n, 3) + • • • 
{z-r-{n- l)/3)5(n, n)))). (2.13) 

Secondly, Use synthetic division to obtain {z)n,-al(z — r), a polynomial of degree < 
n—1, with the remainder S{n,0). Then, evaluate {{z)n,-a/{z—r) — S{n,{)))/{z—r— (3) 
to find the quotient polynomial of degree < n — 2 as well as the remainder S{n, 1). 
Continue this process until a polynomial of degree < 1 left, which is S{n,n—l) + {z — 
r — {n— l)/3)5(n, n). A equivalent description of the above process can be presented 
as follows. Use Homer's method to find 

fir) = (z)„,_„ = 5(n,0) + (z-r)/i(z), deg f,{z) < n - 1, 
where the remainder is S{n,0)). Then, use Horner's method again to evaluate 

fi{z) = 5(n, 1) + (z - r - I3)f2{z), deq f2{z) <d-2, 
which generates the remainder 5(n, 1). Continue the process and finally obtain 

fn-i = S{n, n-l) + {z-r-in- l)/3)5(n, n). 
In short, we obtain S{n,0) = {z)n,-a\^^r' '^('^'1) = {{z)n,-a - S{n,0))/{z - r)|^^^_^^, 

etc. 

Algorithm 2.5 can be demonstrated by the following examples. 
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Example 2.4 For the classical Stirling numbers of the second kind in the case of 
n = 5 and (a, ^,r) = (0, 1,0), from expansion (2.13) we have 

= S{5, 0)+z{S{5, l)+{z-l){S{5, 2)+{z-2){S{5, 3)+{z-3){S{5, A)+{z-A)S{5, 5))))), 
which implies 5(5, 0) =0 and 

z^ = 5(5, 1) + (^ - 1)(5(4, 2) + {z- 2)(5(4, 3) + (^ - 3)(5(4, 4) + (^ - 4)5(5, 5)))). 
Thus, we may use the following division to evaluate 5(5, k) {k = 1, 2, 3, 4, 5). 
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15 
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18 
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25 
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Hence, 5(5, 1) = 1, 5(5, 2) = 15, 5(5, 3) = 25, 5(5, 4) = 10, and 5(5, 5) = 1. 

Prom (2.13) we also immediately know that 5(n, n) = 1 because it is the co- 
efficient of z" on the right-hand side while the coefficient on the left-hand side is 
1. 

Let {ij}j=i be a strictly increasing n-sequence, and let a = {a{j)}j^i be any 
strictly increasing integer sequence in [l,n]. There holds the following well-known 
refinement formula of divided difference (see, for example, [2]) 

(T(fc)-fc) 

f[t,t-t^(i),...,t-t„(^k)]= X! c(j)/[t,ij+i, . . . ,t - ij+fc], 

i=<T(i)-i 

where c(j) = Ct^a > 0. Using this refinement formula one may obtain the refinement 
formula of the generalized Stirling numbers defined by (1.5). 

Proposition 2.6 Let n, /c G No and the parameter triple {a, /3, r) ^ (0, 0, 0) is in 
or . Then there holds refinement formula, 
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^r/(^) ■■= E(-l)' [j] e^'-'^'fi^ + m {z e C), (3.1) 



cr(fc)-fc) 

J=<^(1)-1 

where 

A(A:,)/ := /[i, i + ^/3, t + + 1)/?, . . . , t + j/3] 

3 Generalized Stirling functions 

We now extend the Stirling numbers S{n, k)) expressed by (2.5) to a more wider 
generation form using the idea of [6]. First, in order to cover as large a function 
class as possible, we recall that the generalized fractional difference operator A^'*^ 
with an exponential factor, which is introduced in [6]. More precisely, for ry e C, 
/? € M+, e > 0, the generalized fractional difference operator A^'^ is defined for 
"sufficient good" functions / by 

where (Jl) are the general binomial coefficients given by 

^]=Ml.,= ^(^-'y-(^-^ + ^) (,-eN), (3.2) 
jj .]' .V 

with [/3]o = 1. Noting the generalized Stirling numbers S{n,k) can be represented 
by (2.6), or equivalently, 

which has an extension shown in (2.9). We now extend (2.9) to a more generalized 
form shown in the following definition. 

Definition 3.1 The generalized Stirling functions, ^(7, 77, a, jS, r; e) for any complex 
numbers 7 and rj are given by 

5(7, m e) = S{7, V, a, /3, r; e) := ^.^^^^^^^ 1™ A^'^((^)^,_„) (e > 0), (3.3) 
provided the limit exists; or equivalently, by 



5(7, V, e) = 5(7, V, a, /3, r; e) = (j) ^^'~'^Hr+{v-m^,-a il 7^ 

(3.4) 



provided the series converges absolutely, and 
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From (3.4), we immediately liave 

5(7,0;e) = (r)^,_„ (7 # 0). (3.6) 
Now, an explicit expression of S{-j, r]; e) can be given by the following result. 

Theorem 3.2 // 7 G C and either of the conditions 77 £ C (rj ^ Z), e > 0, or 
r] G C (?7 ^ Z, Re{r]) > Re(^)), e = hold, then the generalized Stirling functions 
5(7, ?7;e) can he represented in the form (3.4) and S{0,r];e) = Sr^fi. In particular, 
ifj — neNo,i] = ke'N, and e > 0, then the corresponding generalized Stirling 
functions S{n,k;e) has the representation (3.4) ■ 

We now present the recurrence relation of the generalized Stirling functions de- 
fined by (3.4) by using the recurrence relations of the generalized raising and falling 
factorial functions shown in Proposition 1.3. 

Theorem 3.3 There hold the following three results. 

(a) For 'J e C, ri e C (rj ^ Z), and e > 0, the generalized Stirling functions 
5(7, 77; e) defined by (3.4) satisfy 

5(7, V, e) = (r + ?7/3 - (7 - l)a)5(7 - 1, V, e) + S{j-l,r]- 1; e). (3.7) 

(b) Let 7 e C, ?7 e C (ry ^ Z), and Re{rj) > Re^y)). The generalized Stirling 
functions S{'y, rf) = S{-y, 77; 0) satisfy 

S{j, 77) = (r + r?/3 - (7 - l)a)S{j + S{j - 1, 77 - 1). (3.8) 

(c) For 7 G C, A: G N, and e > 0, the generalized Stirling functions S{'y,k;e;h) 
defined by (3.4) satisfy 

5(7, k; e) = (r + - (7 - l)a)5(7 - 1, k; e) + S{^-l,k- 1; e). (3.9) 
In particular, 

5(7, k) = {r + kp-{-f- l)a)5(7 - 1, fc) + ^(7 -l,k- 1). 

Clearly, Theorem 6 in [6] is a special case of Theorem 3.3 for a,(3 = 0. And 
Theorem 3, Corollaries 3.1 and 3.2 in [6] are special cases of Theorem 3.3 for a,/3 = 
and 7 = n G N. 

Now we construct the exponential generating function for the generalized Stirling 
functions S{n, 77; e). 

Theorem 3.4 Let z G C, 77 G C, and e > 0. The generating function for the 
generalized Stirling functions 5(7,77; e) defined by (3.4) with 7 = 71 and a/3 ^ is 



^(i+o..r/°(-'"+°;'"°-'V^ 



r{, + 1) 

for 7] and e > 0^ and 



n>0 



(3.10) 
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1 



( 



1 



) 



k 



,n 



(1 + azY^' 



^S{n,k;e) 



z 



(3.11) 



n>0 



for = fc e No and e > 0. 

Remark 3.1 The condition a/3 7^ is not necessary for the left-hand side of (3.10). 
In fact, taking r = 0, /3 = 1, and letting a 0+, we see that (3.10) yields the 
generating function for the generalized Stirling functions of the second kind: 



which was studied in Theorem 4 of [6], and it can be considered as a particular case 
of our Theorem 3.4. 

Similarly, taking e, r = 0, a = 1 and letting ^ — >■ 0+ yields the generating function 
of the generalized Stirling functions of the first kind: 



4 More properties of the generalized Stirling func- 
tions and numbers 

let us consider the set of formal power scries (f.p.s.) J' = 'R[[t;{ck}]] or C[[t; {c}]] 
(where c = (co,ci,C2, . . .) satisfies Co = 1, > for all fc = 1,2, . . .); the order of 
f{t) S J, f{t) = X^feLo fkt'^/^k, is the minimal number r e N such that fr 7^ 0; J'r is 
the set of formal power scries of order r. It is known that J'q is the set of invertible 
f.p.s. and 3^1 is the set of compositionally invertible f.p.s., that is, the f.p.s. 's f{t) for 
which the compositional inverse f{t) exists such that f{f{t)) = /(/(<)) = t. We call 
the element g € J with the form g{x) = X]fc>o ^ generalized power series (GPS) 
associated with {c„} or, simply, a (c)-GPS, and the GPS set associated with {c„}. 
In particular, when c = (1,1,...), the corresponding 3^ and If^ denote the classical 
formal power scries and the classical formal power scries of order r, respectively. 

Wc now develop a kind of asymptotic expansions for the generalized Stirling 
functions S{n, fj,, r; e) = S{n, n, a, l3, r; e) and S{n, fj., jir; e) = S{n, /i, a, /3, /ir; e) and 
generalized Stirhng numbers S{n + ii, fi, r) = S{n+^, /j,, a, (3, r) and S{n + fi, fj,, fir) = 
S{n + fi, fi,a, (3, fir) for large fi and n with the condition n — 0(^^/^) [^l 00). The 
asymptotic expansions of Hsu and Shiue Stirling numbers in [24] and Tsylova Stirling 
numbers in [46], involving a generalization of Moser and Wyman's result [30], are 
included as particular cases. 

The major tool of construction of the asymptotic expansion is the known result 
about the asymptotic formula for the coefficients of power-type generating functions 
involving large parameters shown in [22]. Let a[n) be the set of partition of n 
(n G N), which can be represented by 1*^1 2*^^ . . . j^k^ \k\ -|- 2fc2 nfc„ = n, 

fcj > (j = 1, 2, . . . , n), and with fc = fci + fc2 + • • • + fcn expressing the number of 
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the parts of the partition. For given fc (1 < fc < n), we denote by cr(n, k) the subset 
of cr(n) consisting of partitions of n having k parts. 

Let 4i{z) = X]n>o '^nZ^ be a formal power series over the complex field C in 3^o, 
with qq = g{0) = 1. For every j {0 < j < n) define 

a(n,n—j) 

where the summation is taken over all such partition \^^2^^ ■ ■ ■ n'^" of n that have 
n — j parts. We have the following known result (see for instance [23]): 

For a fixed m G N and for large jj, and n such that n = o{ij}/'^) {jjL — >■ oo), we 
have the asymptotic expansion 

where W{n,j) are given by (4.1). (4.2) is used to derive the Hsu-Shiue Stirling num- 
bers in [24]. We now generalize (4.2) and the corresponding argument to give asymp- 
totic expansion formulas of generalized Stirling functions S{n, /U, r; e) = S{n, /x, a, /3, 
r;e), S{n, n, iJ,r;e) = S{n, n,a, p, fj,r;e), S{n + n,fj,,r) = S{n + fi, fj,,a, /3,r) and 
S{n+ii, n, jjLr) = S{n+iJ,, ji, a, /3, /xr) for large ji and n with the condition n = 0(/ti^/^) 
as /U — >■ cxD. 

Let g{z) = X^„>o o-nZ^ be a formal power series over the complex field C in Jq, 
with ao = g{Q) ^ 0. We may write 

5(.) = aoE7^-"- 

n>0 

For a fixed to <E N and for large /i and n such that n = o{ij}/'^) {jj, — >■ oo), Prom 
formulas (4.1) and (4.2) we have the asymptotic expansion 



where M^(n,j) are given by (4.1). In particular, when n is fixed, the remainder 
estimate becomes 0{^~'^~^). 

To apply (4.2) to the generalized Stirling mmibers 5(7, r}; e) defined by (3.4) with 
^ = n, r] = n and a(3 ^ 0, let us use (3.11) to take 

,(.) = (1 + a.)^/«£!il±i^£)!^^ . J2 ^(!^." (4.4) 

n>0 



when e 7^ 0, and 



n>0 ^ ' 
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when e 0, so that g(0) = (e*^ — (e 7^ 0) and g{0) = 1 (e = 0) not being zero in 
both cases, where S{n, 1; e) = S{n, 1, a, /3, r; e) and S{n + 1,1) = S{n + 1, 1, a, /3, r), 
g{0) = {e^ — Consequently, from (3.11) we have 



{g{z)r = (1 + a^)''''/"^ 



e^(l + a2)'^/" - 1 



/3 



n>0 



for e 7^ 0, and 



(ff(^))'^ = il + az) 



^ (l + az)^/"- l^" 



for e = 0. Therefore, making use of (4.3) yields 

S{n,i^,a,(3,iir; e) 



j=0 



^ J [iJ'-n+jlj [ye^-lj [iJL-n + m]mj 

(4.8) 

for e 7^ 0, and 

S{n + iJ,i^,a,(3,IJ.r) W{n,j) iy(n,w) \ 

[;u]„[n + M]^ + \[n-n + m]mj 

for e = 0, where n = o(/i^/^) as /x — >■ cxd and VF(n, j) (j = 0, 1, 2, . . .) are given by 
(4.1) with Gj being determined by (4.4); namely, for e ^ 0, ao = (e^ — l)//3 and 

aj = [z^]giz) = (4.10) 

while for e = 0, ao = 1 and 

= [zTgiz) = (4.11) 

The coefHcients defined by (4.10) and (4.11) can be evaluated by using the Vandermonde- 
Chu formula as follows. From (4.4), for j = 1, 2, . . . , we have 
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[z^]g{z) = [z^]{l + az) 



r/a 



- 1 ^ f/3/a\ , 



e>o 



e>o k>i 



r/a\ f j3 / a 
k 



(az) 



e+k 



e" - 1 
-1 

1 



r/a 
3 



fe=i 



r/a \ f l3/a 



j — kj \ k 



(r) 



J, -a 



-a^ 



r/a + P/a 
j 



r / a 
j 



= j^[{r + P)j,-a + {e'-2){r)j,.^]. 



Here, the classical Vandermonde-Chu convolution formula we used above, regarded 
as "perhaps the most widely used combinatorial identity" (see P. 8 in [35] by Riordan 
and PP. 51, 61, 64, and 227 in [1] by Andrews), which can be written as 



E 



fe=0 

Similarly, we obtain 



y 

kj \n — k 



x + y 
n 



(a;,y e M,n e No). 



1 



{r)j 



u + m 

for j = 0, 1, 2, Hence, we may survey the above into the following theorem. 

Theorem 4.1 Fore ^ 0, there holds the asymptotic expansion (4-8) ofS{n, fi, fir; e) = 
S{n, II, a, P, /ur; e) for n with n = o{ii^^^) (n oo), where W{n,j) is defined by (4-1) 
with ao = (e^ — l)//3 



«i = ^ [(^ + P)j,-c^ + (e' - 2)(r),-,-„] (j = 1, 2, . . .). 

For e = 0, there holds the asymptotic expansion (4-9) of S{n + fi, fi, fir) = S{n + 
fi, fi,a, (3, fj,r) for n with n = o{fi^/'^) (fi — >■ oo^, where W{n,j) is defined by (4-1) 
with 



{j + m 



[(r + /3)j+i,-c« - {r)j+i-a\ j = 0, 1, . . . . 



Since the formulas (4.8) and (4.9) with W{n,i) and Oj presented in (4.1) and 
Theorem 4.1, respectively, are algebraic analytic identities, we may replace r by r/fi 
in the formulas and obtain the following corollary. 
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Corollary 4.2 For e =^ 0, by replacing the quantity r by r / fi, the asymptotic expan- 
sion (4-8) is also applicahle to S{n, /i, r; e) = S{n, /i, a, (3, r; e) for n with n = o{ij}^'^) 
oo), where W{n,j) is defined by (4-1) with uq = (e^ — 1)//3 and 



1 

1^ 



■13 



+ {e'- - 2) 



0- = l,2,...). 



For e = 0, by replacing the quantity r by r/ pi, the asymptotic expansion (4-9) is also 
applicable to S{n + fi, ji, r) = S{n + fi, ji, a, P, r) for n with n = o{ij}/^) (ji ^ oo), 
where W{n,j) is defined by (4-1) with 



Uj = 



i = o,i, 



In the recent literature, special emphasis has been given to the concept of Riordan 
arrays, which are a generalization of the well-known Pascal triangle. Riordan arrays 
arc infinite, lower triangular matrices defined by the generating function of their 
columns. They form a group, called the Riordan group (see Shapiro et al. [42]). 
Some of the main results on the Riordan group and its application to combinatorial 
sums and identities can be found in Sprugnoli [43, 44], on subgroups of the Riordan 
group in Peart and Woan [33] and Shapiro [39] , on some characterizations of Riordan 
matrices in Rogers [36], Merlini et al. [28] and He et al. [21], and on many interesting 
related results in Chcon ct al. [11, 12], He et al. [19], Nkwanta [32], Shapiro [40, 41], 
and so forth. We now generalize the Riordan arrays associated with classical power 
series to those associated with (c)-GPS, where c = {ck = k\}k>o- The Riordan 
arrays associated with other (c)-GPS can be found in author's later paper. More 
precisely, let c = {cfc = A:!}fc>o. The (c)-Riordan array generated by d{t) £ and 
h{t) G Ji with respect to {ck}k>o is an infinite complex matrix [dn,k]o<k<n, whose 
bivariate generating function has the form 



Fit,x) 



n,k 



k 

n! 



d{t)e 



xh{t) 



(4.12) 



which is called a Shejjer type Riordan array. 

Thus, the (n, k) entry of (c)-Riordan array [dn^k] is 



dn,k 



(4.13) 



for all < A: < n and dn,k = otherwise. It is easy to see that a lower triangular 
array [dn,k] is a (c)-Riordan array if and only if the array {k\dn,k/n\) is a (l)-Riordan 
array, i.e., a classical Riordan array. Evidently, [dn,k] = {d{t), h(t)) can be written 
as 



[d„,fc] = D[[t'^\d{t)(h{t)f\n>k>oD-\ (4.14) 

where D = diag{l, 1, 2!, . . .). 

Rogers [36] introduced tlie concept of the A-sequence for the classical Riordan 
arrays; Merlini et al. [28] introduced the related concept of the Z-sequence and 
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showed that those two concepts, together with the element do.Oi completely charac- 
terize a proper classical Riordan array. In [21], Sprugnoli and the author consider 
the c;liaracterization of Riordan arrays, their multiplications, and their inverses by 
means of the A- and Z-sequences. 

In [36], Rogers states that for every proper Riordan array D = {d{t), h{t)) there 
exists a sequence A = (a;j)fe£N such that for every n. A; e N we have: 

= a^\e\d{t){h{t)f + ai[i"]d(i)(/i(i))'=+i + a2[i"]d(t)(/i(t))'=+2 + • • • 

OO 

= ^a,r]d(i)(/i(i))'=+^- (4.15) 
j=o 

where the sum is actually finite since (i„,fe = 0, VA; > n. We can reformulate it to 

the generalized (c)-Riordan array as follows. 

Theorem 4.3 An infinite lower triangular array D ~ (rfn,fc)„ = {d{t), h(t)) is a 
(c)-Riordan array if and only if a sequence A = (ao ^ 0, ai, a2, . . .) exists such that 
for every n,k gN relation 

«n+l,fe+l — <JO"n,fc H CllUn,k+l H (i2(tn,k+2 + • • ' — djUn^k+j 

Cn+l (^n (^n C^ j— ^ 

(4.16) 

holds. In addition, the generating function A{t) of A— sequence is uniquely deter- 
mined by tA{h{t)) = h{t). 

We now use Theorem 4.3 to establish a new recursive relationship of generalized 
Stirling numbers. Prom expression (3.11) in Theorem 3.4 with e = and a/3 ^ 0, 
we have the generating function of the generalized Stirling numbers shown below: 

l(l + a.)^/"(ii±il^i^)'=|:5(n,A)§. (4.17) 

Theorem 4.4 Let ap ^ 0. The A— sequence (a„)„gN„ of the Riordan array of the 
generalized Stirling number array [dn,k = k\S{n, k) /n\]o<k<n satisfies 

ao = l, a„ = -Lj2^^J^^±l^ (4.18) 

(fc + 1)! 



for all n>l. 



To find the fist column of the array [dn^k]o<k<n, we consider (4.17) for fc = and 
have 

(l + az)'^/" = ^^^^^". 



n>0 



On the other hand, 
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(i+a.r/" = ^(^'/;M(azr 



n>0 - " 

Comparing the right-hand sides of the last two equations, we obtain 

S{n, 0) = S{n, 0, a, /3, r) = n\ (^^^^^ a"" = (r)„,_„. (4.19) 

Formula (4.19) was given in (1.7) and also in (2.9), which are derived by different 
approaches. 

Prom (4.17) we have 



[dn,k] 



0<k<n 



^5(n,fc) 



(4.20) 

0<fe<n 



where S{n,k) = S{n,k,a, (3,r) {aj3 ^ 0). Therefore, surveying the above pro- 
cess, we obtain an algorithm to evaluate generalized Stirling numbers S{n,k) = 
S{n, k, a, /3, r) with a/3 ^ 0. 

Algorithm 4.5 Denote d{t) = {l+azy/" andh{z) = ((l-|-a^)^/"-l)//3 (ajS ^0). 

Let n, fc G No and aji ^ 0. Then we m,ay find A-sequence (a„)„gNo shown in (4-18) 
and establish the array (4-20) except its first column by using the recursive relation 
(4-16) shown in Theorem 4-3, i.e., 

^5(n, k) = Y, - 1, + i - 1) (4.21) 

for all 1 < k < n. The first column of array (4-20) can be constructed by using 
(4-19). Thus, the nth entry of the first column is 

^5(n,0) = -^^^^. (4.22) 
n! n! 

Finally, all S{n, k) = S{n, k,a,p,r) (0 < k < n) can be read from a modification of 
array (4.20); namely from 



— [S{n, fc)]o<fe<„ , 

0<k<n 



where S{n, k) = n J2j>a + J ~ l]j-i'5'('^ ~ 1) ^ + j ~ 1) when 1 < k < n, and 
S{n,0) can be obtained from (4.22) or (4-19). 

Remark 4.1 Similar to the argument in Remark 3.1. the condition afi 7^ in 
Theorem 4.4 and Algorithm 4.5 is not necessary. Algorithm 4.5 can be modified 
to adapt some of cases when a 13 = 0. We will show the application of Algorithm 
4.5 to the calculations of the classical Stirling numbers of the second and the first 
kind, i.e., S(n,k,a, j3,r)) — S{n,k, 0,1,0) and S{n,k,a,j3,r) = S(n,k, 1,0,0), in 
Examples 4.2 and 4.3, respectively. 
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Example 4.1 For the Howard's weighted degenerated Sth'Ung numbers S{n,k) = 
S{n, k, 1, 1, —1). From Algorithm 4.5 or Theorem 4.4, we immediately have gener- 
ating hmetion of the eorresponding A-sequenee A{z) — 1. Then, using (4.21) and 



(4.22) we obtain the Riordan array [d, 



'n,k a<k<r 



[7^S{n,k)] 



0<k<n 



as 



k\ 



S{n, k) 



0<k<n 



Therefore, 



[S{n,k)] 



0<k<n 



which gives S'(0,0) = 1; 5(1,0) 
S'(2,2) = 1; S'(3,0) = -6, 3(3,1) 





1 






-1 1 






1 -1 


1 




-1 1 


-1 1 




1 -1 


1 -1 1 


1 






-1 


1 




2 


-2 


1 


-6 


6 


-3 1 


24 


-24 


12 -4 1 


-1, 


3(1,1) 


= 1; ^(2,0) 



= 6, 5(3,2) = -3, 5(3,3) 



= 2, 5(2,1) = -2, 
1; and 5(4,0) = 24, 



5(4, 1) = -24, 5(4, 2) = 12, 5(4, 3) = -4, and 5(4, 4) = 1 row by row. 

Example 4.2 As we have presented in Remarks 3.1 and 4.1, the condition af3 ^ in 
Theorems 3.4 and 4.4 and Algorithm 4.5 is not necessary. Here, we demonstrate how 
to modify Algorithm 4.5 for the case of (a, 0, r) = (0, 1, 0). The generating function 
of the corresponding classical Stirling numbers {3(n,k) = 5(n, fc, 0, 1, 0)}o<fe<n of 
the second kind is 



1 7" 

-(e-l)'= = E5(n,.)-. 



n>0 



Thus the corresponding Riordan array has generating functions d{z) = 1 and h{z) = 
— 1. Since the compositional inverse of h(z) is h(z) = ln(l + z), the j4-sequence 
characterization of the Riordan array has generating function 



A{z) = 



1 



ln(l 



(-i)'-i fc 

l^k>l k ^ 



E 



(-1)'° 

fe>0 k+1 



which coefficients {a„}„>o, i.e., the elements of A-sequence, can be solved from the 
above equation as 



ao = 1, 



V- (-1)* V- 

"z^ "fcXY = 2^'*" 



n+l 



(-1)* 



fe=l 



k=2 



Thus, we obtain the first few a- 
ao = 1, ai = ^ 



1 1 
02 = - — , as = 



-fe+i- 



a4 = — 



2' 12' 24 

Similar to Algorithm 4.5, we may find the Riordan array 



19 
720' 



(n > 1). 



etc. 
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[dn,k. 



0<k<n = 



S{n,k) 



0<k<n 



S{n, k) 



0<k<n 



1 

1 1 
7 3 



24 12 2 

The Riordan Stirling array of the Stirhng numbers of the second kind is 



[S{n,k)] 



0<k<n 



k\ 



S{n, k) 



0<k<n 



1 

1 

1 1 

13 1 

17 6 1 



which gives all S{n,k) = 5(n,fc, 0,1,0) for < A; < 4. For instance, 5(4,0) = 0, 
5(4, 1) = 1, 5(4, 2) = 7, 5(4, 3) = 6, and 5(4, 4) = 1. 

Example 4.3 For (a,/3,r) = (1,0,0), we can also applied a modification of 
Algorithm 4.5 to evaluate the classical Stirling numbers of the first kind s{n, k) = 
5(n,fc, 1,0,0) as follows. In this case, wc have the corresponding Riordan array 
{d{z),h{z)) = (l,ln(l + z)). Thus the compositional inverse of h{z) = — 1. Thus 
the A-sequence {an}n>o has its generating function 



A{z) = 



1 



Solve the above equation to obtain 

1 



^fe>0 (fc+l) 



ao = 1, ai = — -, 02 = —, as = 0, = —3777: etc. 



12 



720 



which brings us the Riordan array 



[dn, 



k\0<k<n 



k\ 

—s{n,k) 
n! 



0<fe<ra 



fc! 

— s(n,fc) 
n! 



0<;c<ra 



-I 11 

" 4 12 



The Riordan Stirling array of the signed Stirling numbers of the first kind is 



[s{n, k)] 



0<k<n 



-s(n,fc) 



0<k<n 



1 

1 
-1 
2 
-6 



1 

-3 
11 



1 

-6 1 



which gives all s{n,k) = 5(n, fc, 1,0, 0) for < fc < 4. For instance, s(4,0) — 0, 
s(4, 1) = -6, s(4,2) = 11, s(4,3) = -6, and s(4,4) = 1. Of course, the Stirling 
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numbers of the first kind can be evaluated more easily by using formula (2.5) in 
Theorem 2.2, namely, 



s{n,k) ^ S{n,k,l,0,0) = ^ -^[zU 



2=0 



which are simply the coefficients of the powers of z in the expansion of [z]n- 
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